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ABSTRACT  —  Tight  bounds  for  the  Shannon  capacity  of 
amplify-and-forward  cooperative  diversity  networks  are 
derived  for  three  different  adaptive  source  transmission 
policies  in  a  myriad  of  fading  environments:  (i)  constant 
power  with  optimal  rate  adaptation  (ORA);  (ii)  optimal  joint 
power  and  rate  adaptation  (OPRA);  and  (Hi)  fixed  rate  with 
truncated  channel  inversion  (TCIFR).  Our  unified  framework 
based  on  the  moment-generating  function  (MGF)  approach 
allows  us  to  gain  insights  as  to  how  fade  distributions  and 
dissimilar  fading  statistics  across  the  distinct  communication 
links  will  affect  the  Shannon  capacity,  without  imposing  any 
restrictions  on  the  fading  parameters. 

Index  Terms  —  wireless  link  adaptation,  cooperative  diversity, 
channel  capacity 

I.  Introduction 

While  multiple-input-multiple  output  (MIMO)  architectures 
could  drastically  improve  the  range  and  reliability  of  beyond 
line-of-sight  and/or  over-the-horizon  communication  links 
without  increasing  transmit-power  -  the  cost,  weight,  and  poor 
aero-dynamics  of  antenna  arrays  may  prohibit  their  use  on 
both  unmanned  vehicles  and  dismounted  tactical  war-fighters. 
In  multi-vehicle  cooperative  operations,  networked  nodes  in 
a  tight  cluster  may  coordinate  both  their  transmissions  and/or 
receptions  to  mimic  a  space-time  processing  system  as  if  they 
were  part  of  a  single  antenna  array  platform  (e.g.,  [l]-[2]). 
Thus  cooperative  airborne  networking  could  significantly 
increase  the  range  and  reliability  of  a  long-haul  inter-cluster 
communication,  thereby  improving  platform  endurance  with 
enhanced  LPI/LPD  capability,  without  using  an  antenna  array. 
Several  standardization  groups  such  as  IEEE  802.16  and  IEEE 
802.11  have  also  incorporated  cooperative  relaying  into  their 
emerging  wireless  standards  (e.g.,  Mobile  Multihop  Relaying 
Group  has  defined  a  multihop  relay  architecture  in  the  baseline 
IEEE  802. 1 6j  standard). 

Although  an  intermediate  (relay)  node  in  a  cooperative 
wireless  network  may  either  amplify  what  it  receives  (in  case 
of  amplify-and-forward  relaying)  or  digitally  decodes,  and  re¬ 
encodes  the  source  message  (in  case  of  decode-and-forward 
relaying)  before  re-transmitting  it  to  the  destination  node,  we 
shall  focus  on  the  amplify-and-forward  (analog)  relaying 
scheme  in  this  article  because  it  does  not  require 
“sophisticated”  transceivers  at  the  relays.  Nevertheless,  our 
framework  is  also  applicable  for  digital  relaying  once  the 
MGF  of  signal-to-noise  ratio  (SNR)  is  found  or  available. 
While  this  protocol  can  achieve  full  diversity  using  a  virtual 
antenna  array,  there  is  a  loss  of  spectral  efficiency  due  to  its 


inherent  half-duplex  operation.  This  penalty  could  be 
“recovered”  to  some  extent  by  combining  the  cooperative 
diversity  with  a  link  adaptation  mechanism  wherein  the  power 
level,  signal  constellation  size,  coding  rate  or  other 
transmission  parameters  are  adapted  autonomously  in 
response  to  fluctuations  in  the  channel  conditions.  Besides, 
fixed  transmission  methods  that  are  designed  to  provide  the 
required  quality  of  service  in  the  “worst-case”  scenario  are 
very  inefficient  when  better  channel  conditions  prevail. 

But  the  art  of  adaptive  link  layer  in  a  cooperative  wireless 
network  is  still  in  its  infancy  especially  when  optimized  in  a 
cross-layer  design  paradigm.  Majority  of  the  literature  on 
cooperative  diversity  are  limited  to  both  fixed  signaling  rate 
and/or  constant  transmit  power  for  all  communication  nodes. 
For  instance,  while  [3]-[5]  have  studied  the  problem  of 
optimal  power  allocation  in  a  cooperative  wireless  network, 
source  rate -adaptation  [6]  was  not  considered,  and  more 
critically  its  solution  requires  the  knowledge  of  channel  side 
information  (CSI)  of  all  links.  Motivated  by  these 
observations,  [7]  derived  bounds  for  the  Shannon  capacity  of  a 
non-regenerative  link-adaptive  cooperative  diversity  system 
with  limited  CSI,  in  which  the  rate  and/or  power  level  at  the 
source  node  is  adapted  according  to  the  channel  condition 
(i.e.,  only  feedback  of  the  effective  SNR  at  the  destination 
node  is  required  to  be  available  at  the  source  node)  while  the 
relays  simply  amplify  and  forward  the  signals.  However,  their 
analysis  is  limited  to  Rayleigh  fading.  But  in  an  airborne 
platform,  it  is  much  more  reasonable  to  model  the  channel 
gain  of  each  link  as  a  Nakagami-m  or  a  Rice  random  variable 
(due  to  the  increased  likelihood  of  the  presence  of  a  strong 
specular  component). 

The  main  objective  of  this  work  to  develop  a  unified 
framework  based  of  the  MGF  method  for  evaluating  the 
ergodic  capacity  of  cooperative  relaying  systems  with  adaptive 
source  transmission  policies  (since  the  MGF  of  the  end-to-end 
SNR  is  much  easier  to  compute  or  may  be  readily  available). 
This  is  of  significant  interest  because  several  authors  [8]-[10] 
have  recently  argued  that  although  the  MGF  approach  has 
been  successfully  and  extensively  applied  for  evaluating  the 
performance  of  wireless  relaying  systems  in  terms  of  outage 
probability  and  error  rates,  there  have  been  very  limited 
contributions  on  ergodic  capacity  of  fading  relay  channels  [8, 
pp.  2286]  (which  may  be  attributed  to  the  difficulty  in 
evaluating  the  probability  density  function  of  the  end-to-end 
SNR  in  closed-form)  or  explicitly  highlighted  the  complexity 
of  generalizing  the  MGF  approach  for  channel  capacity 
computation  [9].  In  [7],  the  authors  circumvented  this 
difficulty  by  evaluating  upper  and  lower  bounds  on  the 
capacity  instead,  while  [8]  resorted  to  the  method  of  moments. 
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Fig.  1:  Link-adaptive  cooperative  diversity  system  for  ensuring  the  connectivity  and  network  stability  needed  to  support  varying 

quality-of-service  requirements  in  an  airborne  platform. 


Nevertheless,  both  their  analyses  and  results  are  limited  to 
Rayleigh  channels.  In  [10],  the  authors  proposed  a  general 
method  for  channel  capacity  analysis  using  an  integral  relation 
known  as  Ej-transform,  but  their  solution  also  requires  the 
derivatives  of  the  MGF  and/or  its  auxiliary  function,  which 
can  be  cumbersome.  Moreover,  the  capacity  formula  for  the 
OPRA  policy  is  an  approximation  while  the  capacity  for  the 
TC1FR  policy  was  not  considered.  More  recently,  we  have 
presented  yet  another  unified  approach  for  ergodic  capacity 
computation  based  on  a  cumulative  distribution  function 
(CDF)  method  [11]  (which  utilizes  the  MGF  of  end-to-end 
SNR  in  conjunction  with  a  fixed-Talbot  method  [12]). 

II.  System  Model 

In  Fig.  1,  the  source  node  S  communicates  with  destination 
node  D  via  direct-link  and  through  N  amplify-and-forward 
relays,  R;/e  {1,2 in  two  transmission  phases.  During 
the  initial  Phase  I,  S  broadcasts  the  signal  x  to  D  as  well  as  to 
the  relays  R ,  where  the  fading  coefficients  between  S  and  D,  S 
and  the  /-th  relay  node  //„  and  R,  and  D  are  denoted  by  a  d , 

as ..and or. d ,  respectively.  The  received  signals  at  the 
destination  and  the  relay  is  denoted  respectively  by 


y,A 

=  JP^as,aX  +  risd 

(1) 

y,4 

=  Jp~  ax  +  T)s. 

(2) 

where  Psd  =  PSI  is  the  transmitted  power  at  source,  while  i)s  d , 
>lsi  are  the  additive  noises  introduced  between  source  to 

destination,  and  source  to  relay,  respectively.  During  the 
second  phase  of  cooperation,  each  of  the  N  relays  transmits  the 
received  signal  after  amplification  via  orthogonal 
transmissions  (e.g.,  TDMA  in  a  round-robin  fashion  and  /or 
FDMA).  Suppose  Pid  is  the  transmit  power  of  the  /-th  relay 

while  transmitting  its  signal  to  the  destination  node.  The 
received  signal  at  the  destination  during  the  second  phase  of 
operation  is  given  by  [1] 

y,4 =  GiCci,dySj +r/ij  (3) 

where  G,  is  the  /-th  relay  amplifier  gain,  denoted  as  [1], 

Ch  =  Jp~/fp~\  as4\~  +N0  ,  T)id  is  the  noise  introduced 


between  /-th  relay  and  destination.  We  further  assume  that  the 
total  transmitted  power  Psd  +  f^Pd  =  Pr  is  fixed. 

i  =  1 


Suppose  if  a  maximum-ratio  combining  (MRC)  is 
employed  at  the  destination  node,  the  total  received  SNR  can 
be  shown  as  [1] 


/v 

7t  =  7s4  +  Yj 


r.tYiA 


1+r.;  +  r,A 

where  ysi  =  |«.  PsIt  /  N0  and  yia 


TV 

■■7,a+YJi  =  7t*  (4) 


'  Pi  drs  i  N0  denote 


the  instantaneous  SNRs  of  the  respective  links,  while  the 
instantaneous  SNR  of  a  two-hops  path  can  be  accurately 
approximated  as  the  harmonic  mean  of  the  link  SNRs.  The 
total  SNR  may  be  upper  and  lower  bounded  as 


where  y  =  min(ysi,y,d) . 

If  Ys,d’Ys,i’Yi,d  are  independent  random  variates,  then  it  is 
straight  forward  to  show  that  the  MGF  of  yrB  is  given  by 

(6) 

/  =  i 


For  instance,  the  MGF  of  Y,  =  YslYtdl(Ysl  +  Y,4)  for  a 


Rayleigh  channel  with  independent  but  non-identically 
distributed  (i.n.d)  fading  statistics  is  well-known,  and  it  is 
given  by  [14] 

(f>7_  (S)  =  [(i  /  n  -  i  /  q  „  )2  +  (i  /  a  .  +  i  /  aid  >]  /  a2 


2s 


A3n  q 


-In 
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5  +  A  +  - 


12  Q 


£2  Q 


(7) 


where  £2ab  =  E[yab\  corresponds  to  the  mean  link  SNR,  and 


A  =  ^(1  /  £2~  - 1  /  £1 ~  )2  +  2s  (1  /  Q , ~ + 1  / 12,  „ )  +  s : 2 


Moreover,  the  MGF  of  yi  for  a  Nakagami  channel  with 
independent  identically  distributed  (i.i.d)  fading  statistics  is 
given  by  [15] 

(  I  Qs) 

fZt(s)=  1FAm,2m\m  +  —  --—\  (8) 

'  V  2  4;w ) 
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If  a  closed-form  expression  for  y.  is  not  available  (e.g., 
Nakagami-m  channel  with  i.n.d  fading  statistics),  but  does 
exist  for  y  =  min  ( Ysf,  Yu) ,  we  may  then  develop  capacity 
bounds  using  the  inequality  used  in  [7],  (i.e.,  see  (5))  viz., 

N  N 

(s)<<Py(s)<<PrJs)Y[<Pfi(sl2)  (9) 

1  =  1  i  =  l 

In  this  case,  the  MGF  of  y.  for  2-hops  path  may  be  derived  as 


4>fis)=  ^  [  e  a  fyk  (*)[1  -  F  (*)]  dx 


(10) 


k  e  {( sj),(i,d )} 
j*k 


For  instance,  it  is  not  difficult  to  show  that  the  MGF  of  y.  for 
Nakagami-m  fading  with  i.n.d  fading  statistics  is  given  by 


f(mk  +  nij  ) 


^  *  <=  m),M)  r(mk)r(mj ) 

j*k 


Q  jmk 


sQ.pk  +Q.jmk  +  Q.kmj  J 


(11) 


x  -^2F[ 
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l—m,mk;  1  +  mk; 


(sQ.k  +mk)Q.j 


sQ.jQ.k  +  Q.jmk+Q.kmj  J 


Once  the  MGF  of  y  or  y.  is  found,  we  can  compute  the 
outage  probability  (i.e.,  it’s  CDF)  efficiently  using  a  fixed 
Talbot  (i.e.,  multi-precision  Laplace  transform  inversion) 
method  [12],  viz., 


+——<!>  0')e", 
2  M 


1  +  j°(.6k), 

■m) 


(12) 


where  r  =  2M /(5x),  a(6k)  =  0t  +  ( 0t  cot( 6t ) - 1 ) cot(<9, ), 


Qk  =  knl M  and  s(8k)  =  rdk (j  +  cot( 6t )). 

In  Section  III,  we  will  show  that  the  ergodic  capacity  of 
ORA,  TCIFR,  and  OPRA  source  adaptive  transmission 
policies  can  be  expressed  in  terms  of  the  MGF  alone.  Hence 
the  ergodic  capacity  may  be  evaluated  readily  for  all  of  the 
cases  discussed  above  (since  the  MGF  of  y  or  y.  is  available 
in  closed-form).  For  instance,  (15),  (17)  and  (22)  in 
conjunction  with  (6),  (9),  (11)  and  (12)  generalize  the  results 
in  [7]  to  i.n.d  Nakagami-m  channels.  Similarly,  precise 
estimates  of  the  ergodic  capacities  with  different  adaptive 
source  transmission  techniques  in  i.n.d  Rayleigh  fading  [or 
i.i.d  Nakagami-m  fading]  can  be  obtained  by  substituting  (6) 
and  (7)  [or  (8)]  into  (15),  (17)  and  (22). 


III.  Ergodic  Capacity  Computation  in  F ading  Channels 

The  well-known  Shannon-Hartley  law  tells  us  that  there  is 
an  absolute  limit  on  the  error-free  bit  rate  R  that  can  be 
transmitted  within  a  certain  channel  bandwidth  B  at  a  specified 
SNR.  This  theoretical  limit  denotes  the  channel  capacity  C. 
Shannon’s  noisy  channel  coding  theorem  also  states  that,  it  is 
not  possible  to  make  the  probability  of  error  tend  to  zero  if  R  > 
C  with  any  code  design.  Thus,  it  is  clear  that  the  metric  C 
plays  an  important  role  in  the  design  or  appraisal  of  any 
communications  system  (since  it  serves  as  an  upper  limit  on 
the  transmission  rate  for  reliable  communications  over  a  noisy 
communication  channel).  In  this  section,  we  will  derive 


generic  expressions  for  computing  the  ergodic  capacity  with 
different  adaptive  source  transmission  techniques. 


A.  ORA  Policy 

When  only  the  rate  is  adapted  by  changing  channel 
conditions,  the  channel  capacity  is  given  by  [6] 


r 

n\ 


f”ln(i+ y)fr(y)dy  (13) 

A  + 1  In  2  J  o 


-'ORA 

B  N  + 1 

Utilizing  the  “exponential  type”  integral  representation  of 
ln(/+l)  (see  Appendix  A),  we  can  facilitate  the  averaging 

problem  in  (13)  given  that  the  MGF  </>y(-)  is  available  in  closed 
form.  Substituting  (A. 5)  into  (13),  we  obtain 


Cn 


B 


I _ I 

A  +  l  In  2 

_ 1 _ 1_ 

N  + 1  In  2 


(i  ^-2'r)fr(r)dr 


rqr 

Jo  X  |_Jo 

J  - — [l-0r(2x)]<fe 


dx 


(14) 


Substituting  y  =  2x,  dy  =  2dx ,  (14)  can  be  re-stated  as 


C, 


1  1 


(15) 

B  N  + 1  In  2  J  “  y  L  ',""J  ' 

Note  that  ^  (•)  is  a  monotonically  decreasing  function  with 
respect  to  its  argument  and  0  <  0  {s)  <  1 . 


B.  TCIFR  Policy 

In  CIFR  policy,  the  transmitter  adapts  its  power  to  maintain 
a  constant  SNR  at  the  receiver  and  uses  fixed-rate  modulation 
and  fixed-code  designs.  This  technique  is  the  least  complex  to 
implement  given  that  reliable  channel  estimates  are  available 
at  the  transmitter.  However,  when  the  channel  experiences 
deep  fades,  the  penalty  in  transmit  power  requirement  with  the 
CIFR  policy  will  be  enormous  because  channel  inversion 
needs  to  compensate  for  the  deep  fades.  To  overcome  this,  a 
truncated  channel  inversion  and  fixed  rate  policy  (TCIFR)  was 
considered  in  [6]  where  the  channel  fading  is  only  inverted 
above  a  fixed  cutoff  fade  depth  y0.  The  data  transmission  is 
ceased,  if  jHalls  below  y0.  In  this  case,  it  is  easy  to  show  that 


the  outage  probability  is  Pout  =Fy(y0)  and  the  channel  capacity 
is  given  by  [6] 


CL 


N  + 


-log. 
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r'fyWr 

'>0  J 


Fy(ylt)  (16) 


-1 1  l 

Substituting  fy(y)  =  —  J  R Q\j>r(-jco)e~m^dco  (i.e.,  i 


inverse 


Fourier  Transform  of  the  characteristic  function)  in  (16),  we 
obtain 


C„ 


1 


B 


N  + 1 


log2|  1-^  I FUy0) 


(17) 


where  V  =  —  f  R e{</>  (-joj)Ei(-jcoy0)}d(o  while  exponential 

7TJo 

r-  e~Ia 

integral  Ei{-jc)  =  J  - dt  can  be  evaluated  using  MATLAB 

as  cosint(c)-  y(#72  +  sinint(c)).  Observe  that,  if  the  MGF  of 
y  is  known  in  closed-form,  then  the  above  integral  can  be 
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evaluated  efficiently  via  Gauss-Chebyshev  quadrature  method 
over  a  wide-range  of  fading  channel  models  and  diversity 
combining  techniques  employed. 


C.  OPRA  Policy 

In  OPRA  scheme,  the  transmission  power  and  rate  is 
matched  to  the  varying  channel  condition  through  use  of  a 
multiplexed  multiple  codebook  design.  This  leads  to  the 
highest  achievable  capacity  with  CSI.  From  [6],  we  have 


C.„ 


——I 

f  +  1  ln2  J  r. 


In 


r_ 

r«j 


fXY)dy 


(18) 


B  N- 

Substituting  (A.4)  into  (18),  and  by  rearranging  the  integral 
terms  we  obtain. 


1  r°°  1  T  r°°  r°°  ~2xy> 

T— [  -\e~2x\  fy(r)dy-  \  e  n  f  (y)dy  dx  (19) 
1  m2  Jo  x  L  J/o  Jn> 


C  1 

'-'OPRA  1 

B  N  + 

which  can  be  re-stated  as 


C, 


OPRA 

B 


since 


=-^Tli-["-[e-2l[i-^(n)]-^(2x/j'o,ro)>  (2°) 

N  +  \  In 2  Jo  x L  7  7  J 

f  fr(Y)dy  =  l-Fr(yo)  =  Fy(yo)  is  the  complementary 

J  Yo 


-2  xy 


CDF  of  y  and  [  e  n  fr(y)dy  =  </>r(2 x/y0,y0)  denotes  the 

J  Yo 

marginal  MGF.  Since  the  desired  marginal  MGF  is  typical  not 
available  in  closed-form,  we  circumvent  this  difficulty  by 
computing  this  term  via  the  CDF  of  an  auxiliary  function 
(details  shown  in  Appendix  B),  viz., 

e~lxv-W\~t[— ]+W  *  (21) 

B  N  +  l  ln2  Jo  x  \y„J 

In  (21),  </>r(2x/yo)  is  available  in  the  closed  form  while  the 


CDF  terms  F  (yj  and  F(y.)  are  compute  efficiently  using  the 

fixed-Talbot  method  (see  (12)). 

To  achieve  the  capacity  shown  in  (21),  the  channel  fade 
level  (i.e.,  CSI)  tracked  at  the  receiver  must  be  conveyed  to 
the  transmitter  on  the  feedback  path  for  power  and  rate 
adaptation  in  accordance  with  the  prevailing  channel 
condition.  When  y<y0,  no  data  is  transmitted,  and  thus  the 
optimal  policy  suffers  an  outage  probability  given  by 
Pm,=Fr(y0),  which  equals  to  the  probability  of  no 
transmission.  The  optimal  cutoff  SNR  must  satisfy 

f;  (To)- r„[i + f  y'fr{y)dy\  =  o  (22) 

J  ro 

The  integral  term  in  (22)  can  be  evaluated  efficiently  similar 
to  the  development  of  (17)  when  the  MGF  is  available. 
Furthermore,  asymptotic  analysis  of  (22)  shows  that  y0  —>  0 

when  the  mean  SNR  Q.  — »  0  because  F(x)  — » 1  and  f(x)  — »  0 
(the  effect  of  Q  — >  0  can  be  predicted  from  the  normalized 
PDF  or  the  normalized  CDF  curve  when  its  argument  x  — >  °°  ). 
When  LI  — >  °°  ,  F  (x)  —>  0  because  this  is  equivalent  to 
computing  the  CDF  when  its  argument  x  — >  0 .  It  is  also  well- 
known  that  0  (•)—»()  as  £2— ><x>.  Hence,  yo  — >  1  as  Q— 
Thus  yo  (can  be  determined  by  solving  (22)  numerically) 
always  lies  in  the  interval  [0,  1]  regardless  of  the  fading 
channel  model  and  number  of  relay  nodes  employed. 


IV.  Computational  Results 

In  this  section,  selected  numerical  results  are  presented  for 
the  Shannon  capacities  of  cooperative  analog  relaying  systems 
under  different  source  adaptive  transmission  policies.  The 
following  parameters  (arbitrarily  chosen)  will  be  used  to 
generate  the  plots,  unless  stated  otherwise:  Qs  l  =  E/N0,  £2  .  = 

0.5 EM),  =  0.5 E Mo,  (-Kj  =  E*/No,  and  Q  d  =  0.2E/N0. 

Fig.  2  depicts  the  comparison  channel  capacities  of  three 
distinct  source  adaptive  transmission  policies  in  an  i.n.d 
Rayleigh  fading  channels  with  two  cooperative  relays.  As 
anticipated,  there  is  no  significant  difference  observed  in  the 
capacities  of  OPRA  and  ORA  at  high  SNR.  The  capacity  of 
TCIFR  is  plotted  for  the  cut-off  SNR  y0=  6  dB.  Although  not 
shown  in  this  figure,  we  also  noticed  that  the  curves 
corresponding  to  the  tight-bound  case  are  in  good  agreement 
with  the  Monte  Carlo  simulation  results.  Moreover,  the  actual 
capacity  is  slightly  closer  to  the  lower  bound  (rather  than  the 
upper  bound)  at  low  SNRs.  Although  the  authors  in  [7]  have 
studied  the  channel  capacities  of  cooperative  relaying  system 
in  i.n.d  Rayleigh  fading  channel,  but  their  framework  does 
lend  itself  to  the  analysis  of  the  “tight-bound”  case  or 
generalize  to  other  fading  channels,  whereas  our  framework 
encapsulates  all  these  cases  in  an  unified  way. 


Fig.  2  Channel  capacities  of  ORA,  OPRA  and  TCIFR  in  an  i.n.d  Rayleigh 
fading  channel  with  two  cooperative  relays. 

Fig  3  illustrates  the  TCIFR  channel  capacities  for  E/N0  = 
6dB  and  EMo  ~  15dB  as  a  function  of  cut-off  SNR  yn  in  an 
i.n.d  Rayleigh  fading  channel  with  two  relays.  It  is  evident  that 
there  exists  an  optimal  choice  of  the  cut-off  SNR  which  will 
maximize  the  channel  capacity  while  EMo  is  fixed. 

Fig.4  shows  the  channel  capacities  (i.e.,  the  tight-bound 
case)  for  different  source  adaptive  transmission  schemes  over 
i.i.d  Nakagami-m  channels  (fading  severity  index  m  =  0.5,  1, 
1.5  and  2).  It  is  apparent  that  the  channel  capacity  increases 
with  the  increasing  value  of  m  (i.e.,  as  channel  experiences 
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less  severe  fading).  Moreover,  the  gap  between  the  ORA  and 
OPRA  capacities  widens  as  m  decreases. 


Fig.  3  TCIFR  channel  capacity  versus  cut-off  SNR  in  an  i.n.d  Rayleigh  fading 
channel  with  two  relays. 


Fig.  4  Channel  capacities  of  ORA,  OPRA  and  TCIFR  policies  in  an  i.i.d 
Nakagami-m  channel  with  a  single  cooperative  relay. 

Fig.  5  depicts  the  ORA  channel  capacity  for  i.n.d 
Nakagami-m  channels  with  varying  number  of  relays  (N  =  0, 
1,  2)  and  transmit  power  allocation.  Without  any  loss  of 
generality,  we  can  show  that  the  average  SNRs  on  different 
links  are  given  as  follows: 


a ,  = 


f  d.  v" 


V  d*4  J 


r  ,  \~" 

“l  4 


d  . 

V  s-d 


cSslET/No,£lsa  = 
cSuET/N,Q2J  = 


d 


V" 


\d'4j 

r  r  ^ 
dij 


\d’4j 


cSmJIEt/N', 

n 

cS2,Et/No 


f dsd A 

and  -  -j- 

{ds,j 


c  Ss  d  Et/No,  where  c  is  constant  that  i 


is 


related  to  the  carrier  wavelength  (i.e.,  c  =  -  for  free 


A 

4n 


space  path  loss),  n  is  the  path  loss  exponent,  Ss  d  =  8  .  =  —L 

PT 

while  8id  =  (1  -  8SJ )  =  since  ^  8a)t  =  1 ,  and  dab 
rT 

denotes  the  distance  between  link  a-b.  To  generate  Fig.  5,  we 
have  arbitrarily  chosen  dsl  =ds2  =dld  =d^d  =500 m ,  c  =  10" 

2,  dsd  =  1000  m ,  and  n  =  4. 

It  is  evident  that  equal  power  allocation  strategy  may  be 
reasonable  but  not  optimal  in  all  cases  (i.e.,  optimum  power 
allocation  strategy  is  more  beneficial  when  the  channel 
experience  more  severe  fading).  It  is  also  observed  that  for  the 
chosen  set  of  parameters,  increasing  the  number  of 
cooperating  relays  beyond  a  certain  limit  (in  this  case,  N  =  1) 
does  not  necessarily  translates  into  higher  spectral  efficiency 
(although  diversity  gain  may  be  harnessed)  -  there  exists  an 
optimum  number  of  cooperating  nodes.  For  instance,  the 
single  relay  case  yields  the  highest  capacity,  followed  by  the 
two  relays  case,  and  finally  the  no  relay  case  for  values  of  8si 

below  0.9  when  m  =  3  and  E7/N0  =  20dB.  We  also  observe 
that  the  tightness  of  upper  and  lower  bounds  on  the  Shannon 
capacity  becomes  looser  as  m  increases.  By  comparing  the 
curves  corresponding  to  the  single  relay  and  no  relay  cases,  we 
can  deduce  that  the  optimum  power  allocation  is  not  strongly 
influenced  by  the  source-destination  link  but  rather  on  the  link 
qualities  between  source-relay  and  relay-destination. 


Fig.  5  Shannon  capacity  of  a  cooperative  relaying  system  consisting  of  N 
relays  with  ORA  policy  but  with  different  fixed  transmit  power  allocations 
in  a  i.n.d  Nakagami-m  channels. 
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Appendix  A 


In  this  appendix,  we  provide  a  sketch  of  our  derivation  for 
an  “exponential-type”  integral  representation  for  In  y .  Such  a 
representation  will  facilitate  the  averaging  problem  that  arises 
in  the  channel  capacity  analysis,  and  therefore  leads  to  a 
unified  approach  for  evaluating  the  ergodic  Shannon  capacity 
in  a  myriad  of  fading  environments  and  relaying  strategies. 

Utilizing  [13,  eq.  (1.512.2)],  we  have 

ifr ziT 


\ny=2Yj- 


■2YJ\{y)k’  r>o  (A.i) 


k  =  l 
k  odd 


where  y  = 


r- 1 


Substituting  yk  = : 


y+ 1  “  -  T(k) 


-f 

(k)  Jo 


xtlellydx  [13, 


eq.  (3.381.4)]  into  (A.I),  we  obtain 


In  y=  lY  -  — —  xt  le  xlf  dx  =  2  f~e  x/y 

4rik  TO)  Jo  Jo 


k  =  1 
k  odd 


k  i  r(k) 


L-i 

k  =  1  H  • 


dx 


(A.2) 


1  —  g  A  y  ^  | 

Recognizing  that  —  s h  x  = - =  >  — xk  [13,  Eq. 

r  v  r—J  lr  I 


(1.411.2)],  (A.2)  can  be  re-stated  as 

In  y  =  2  [  e*ty  —  ( ex  -  e  x  )e 

Jo  x 


[r-l>dx 


(A3) 


dx 


Finally  using  variable  substitution  x  =  z(y—  1),  dz  = — - — ,  we 

y- 1 

arrive  at  (A.4)  after  some  routine  algebraic  manipulations: 

\ny=[  -\  elz -e1:r^dz,  y>  0  (A.4) 

Jo  z  L  J 

It  is  also  obvious  that 

ln(y+Y)—["— — I”  l  —  e~2ir~\dz,  y>- 1  (A.5) 

jo  7  L  J 


Appendix  B 

Let  (j),: (s)  =  [  e  "  fx(x)dx  and  O  v (jo)  =  j  ejmfx(x)dx 

Jo  Jo 

denote  the  MGF  and  the  characteristic  function  (CHF)  of 
random  variable  X  >  0  respectively.  Thus  the  CHF  is  related 
to  the  MGF  as  O  ,  (jo)  =  (j>x  (-  jo)  . 

Also  while  evaluating  the  ergodic  capacity  for  OPRA 
policy,  we  must  compute  the  marginal  MGF  of  SNR.  Suppose 
the  MGF  0  (■)  is  known  in  closed,  we  may  then  use  Abate’s 

fixed  Talbot  method  (multi-precision  Laplace  inversion)  for 
computing  the  marginal  MGF  efficiently.  A  sketch  of  this 
derivation  is  highlighted  next: 

Let  us  define  an  auxiliary  function  fjx)  =  e  /hfr(x) .  Hence, 
(pr(Jj, a)  =  J  e», fr(y)dy=  F.(oo)- F.(a)  (B.l) 

where  FAy)  =  [  fJx)dx  =  [  e“’f  ( x)dx  .  It  is  obvious  that 

FM  =  [el<‘fr(x)dx  =  m  (B.2) 

Furthermore,  we  have 


=  J  e  '\f.(x)dx, 

=  ff(x)dx,  (B.3) 

=  <pr(s  +  p) 

Substituting  (B.2)  and  (B.3)  (in  conjunction  with  (12))  into 
(B.l),  we  arrive  at  a  numerically  efficient  method  for 
evaluating  the  required  marginal  MGF. 
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